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Abstract 



Standard decoupling of heavy fermions may fail when there are non-perturbative 
variations in a scalar field which gives masses to the fermions. One situation of 
phenomenological relevance is the case of sphalerons in the presence of fermions 
at finite temperatures. The free energy of a simple model is determined using 
a non-perturbative technique to study the effect of fermions on the scalar field. 
The effects of quantum and thermal fermionic fluctuations on the free energy of 
the thermal sphaleron are calculated, including contributions from the gradients 
of the scalar field to all orders. 
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1 Introduction 



When a particle is made heavy by increasing a dimensionful parameter, its virtual effects on 
physical quantities are suppressed by inverse powers of the parameter. Under these condi- 
tions the particle is said to decouple QXQ . However, decoupling fails when a particle is made 
heavy by increasing a dimensionless coupling. This is the case for instance when a fermion 
mass is made heavy by increasing a Yukawa coupling, while the Higgs vacuum expectation 
value is held fixed. Therefore, chiral gauge theories in particular can be sensitive to the 
short distance physics associated with heavy fermions, as a consequence of this perturbative 
non-decoupling [0] . 

Recently, Banks and Dabholkar have pointed out that decoupling may also fail when 
there are non-perturbative variations in a scalar field which gives masses to fermions ||. More 
specifically, these authors noted that the effects of heavy fermions cannot be summarized by 
a local effective Lagrangian for the light degrees of freedom if the scalar field which gives 
the fermion its mass goes through a zero. Heuristically, this non-perturbative non-decoupling 
can be understood as stemming from the local vanishing of the fermion mass. To see this 
more clearly, consider the fermion determinant in a scalar background 4>(x). The effective 
action for the the scalar field may be expressed as 

S eff [<f>] = S[4>] - i\ndet[ifl - g( f>] , (1) 

where S [0] is the classical scalar field action, i p is a gauge covariant derivative and g is a 
Yukawa coupling. To evaluate the determinant in an arbitrary background one may attempt 
a gradient expansion. However, any such approximation is surely doomed. The leading order 
term in such an expansion is given by {d^(j)) 2 / g 2 <p 2 - It is clear that this approach fails in the 
case of a vanishing scalar field. Thus, the effects of a heavy fermion can only be summarized 
by a non-local effective lagrangian. 

This non-decoupling can be recognized by an apparent paradox |3|] . Begin by considering 
a chiral gauge theory with heavy mirror fermions for all light fermions. The fermion number 
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currents for the light fermions and for the mirror fermions are assumed to be anomalous, while 
their difference is exactly conserved. Furthermore, assume that at energies small compared 
to the heavy fermion mass, there exists a local effective lagrangian which describes the 
physics of the light degrees of freedom. Now consider the decay of a light fermion through 
the fermion number anomaly. In the constrained instanton background M, the scalar will 
necessarily pass through zero, and furthermore, all of the flavor eigenstates will have zero 
mode solutions to the Dirac operator. In the absence of mixing between heavy and light 
fermions, a non-vanishing decay amplitude necessarily requires heavy fermions on external 
legs to soak up their zero modes. But, this contradicts the initial assumption that a local 
effective theory incorporates all of the physics of the heavy fermions. 

This paper will consider the effects of heavy fermions on background scalar fields with 
non-perturbative variations. Specifically, the effect of heavy fermions on sphalerons at finite 
temperature will be addressed 

It is by now widely accepted that fermion number violation occurs at finite temperatures 
due to the production and decay of sphalerons |6|, [F]. Sphaleron processes may then play 
an instrumental role in determining the baryon asymmetry in the universe fg]. Considering 
the discussion given above, it is not hard to imagine that heavy fermions may influence the 
production and decay of sphalerons. In the process of sphaleron production and decay, the 
fermion number of each flavor of fermion must be violated, as a consequence of the fermion 
number anomaly discussed in the case of instantons above. These fermions may receive 
masses exceeding the mass of the classical sphaleron through Yukawa couplings to the scalar 
field, since the sphaleron mass is determined at tree-level from the bosonic sector of the theory 
only. If the heavy fermions do not mix with lighter fermions, then the sphaleron process 
from a state with no heavy fermions to a state with heavy fermions is prohibited strictly 
on the basis of kinematics^ Setting aside for the moment the obvious technical problems 

^dn realistic models with fermionic mixing, the rate of sphaleron decay will be suppressed by a Cabbibo 
angle which could be quite small for heavy fermions. 
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accompanying a large Yukawa coupling, this simple scenario illustrates that fermions with 
large Yukawa couplings to scalars are expected to modify the scalar background field 0. 

One can imagine at least two ways in which the scalar field could behave under the 
influence of heavy fermions. It may be that the sphaleron is stabilized by the effects of 
heavy fermions, if its negative modes are turned positive by fermion loop effects Q. It may 
also be that the sphaleron mass is simply boosted by effects of heavy fermions, but remains 
an unstable saddlepoint. This latter effect will be demonstrated in the following pages. 

The sphaleron rate in thermal equilibrium, T, is controlled by the Boltzmann exponential 

1 

T oc C x exp[-(3E sph \ (2) 

for temperatures T = 1/(3 below the mass of the sphaleron, E sp h. The sphaleron is a solution 
of the time-independent classical equations of motion for scalar and gauge fields determined 
by minimizing the energy functional. In the one-loop approximation, the prefactor C is 
determined by calculating the gaussian and zero-mode fluctuations around the sphaleron. 
This factor represents the entropy associated with the sphaleron, and it combines with the 
energy to give the free energy of the sphaleron. However, this free energy will not in general 
be a minimum of the free energy functional. The self-consistent method for determining 
the sphaleron rate is to first calculate the free energy functional and then minimize it. This 
determines a new configuration, the thermal sphaleron. The rate is then given by the equation 

T oc exp [-(3 F sph ] . (3) 

This rate will include all of the effects of fermionic fluctuations. 

In this paper such a self-consistent procedure is performed in a simple model, in the limits 



where the fermion mass is greater than or less than the temperature p"T|, |T2"fl . The fermionic 



determinant is found for an arbitrary scalar background, using a WKB approximation, so 



2 A situation like this has been considered recently in an investigation of sphaleron stability in the presence 
of zero- modes Mm. 
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that the effects of gradients to all orders are included. Thus, the problems involved in 
non-locality mentioned above may be addressed. 

In the high T limit, it is shown that the finite temperature effects serve to make the 
gradient expansion well-defined, as expected. This allows a comparison of the WKB result 
with the gradient expansion. The fermion mass is effectively controlled by the temperature, 
and the gradients are suppressed by inverse powers of this dimensionful quantity. However, 
the fermions do not decouple in the high T limit, since they contribute to temperature- 
dependent renormalizations of the parameters in the theory which are physically observ- 
able. In the simple model considered here, the only effect of the fermions in this limit is a 
temperature-dependent renormalization of the scalar mass. In the low T limit, it is shown 
that the heavy fermions also do not decouple, however their effects cannot be summarized by 
a temperature-dependent renormalization of the bare parameters. Furthermore, it is shown 
that the gradients are still suppressed in the regime where quantum fluctuations become 
important. 

The rest of the paper is structured as follows. Section 2 contains a review of the 1+1 
Abelian Higgs model which has many of the features of the electroweak theory necessary for 
a study of finite temperature fermion number violation. In Section 3, the finite temperature 
effective action for the 1+1 Abelian Higgs model is calculated to all orders in gradients, 
using a WKB approximation. The details of these calculations are lef tfor the appendix. 
Section 4 contains a discussion of the contribution of heavy fermions to the free energy, and 
calculations of the thermal sphaleron in high and low temperature regimes. Finally, the last 
section contains a discussion of the results. 

2 The Model 

The 1+1 dimensional Abelian Higgs Model, chirally coupled to fermions, contains several 



features necessary for a study of finite temperature fermion number violation [13, 14, 15 . 
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The Minkowski action of the model is given by: 

S = N J d 2 x ~\fI v + l^l 2 - A (|$| 2 - v 2 /2) (4) 
N f 

+ x: / d 2 xr L iD L r L + r R iD L r R + g$wk + g* *t>i ■ 

a J 

The U(l) gauge field is chirally coupled to A" flavors of fermions through the covariant 
derivative, Dl,r = 7^(9^ ± ieA^) . The complex scalar $ is responsible for spontaneously 
breaking the U(l) symmetry. The dimensionful parameters in the theory are A with dimen- 
sions of mass-squared, and e and g with dimensions of mass. The scalar mass is = \/2\v 
and the fermion mass is m/ = gv 0. By a rescaling of the scalar field, it can be shown that 
the limit of weak scalar self-coupling corresponds to v » 1. 

The work presented in the following sections is based on a large N expansion of the free 
energy of this model so that the effects of a large Yukawa coupling can be considered P, |l2j. 
An overall A" has been included in the bosonic part of eq.^ to make this expansion well- 
defined 0. The large A" limit is taken with the parameters e, g, A and v in eq.|4] held fixed. 
Then, the tree-level bosonic action and the one-loop fermion determinant are the leading- 
order O(N) contributions to the free energy. Bosonic loops are suppressed by powers of 1/N 
and are neglected. 

The model contains a conserved gauge current, J M = yj a i^lI^l ~ ^rIi^r-i an d global 
fermion number currents associated with each fermion flavor, J" = ^£7//0£+V>#7^'0|j, which 
obey the anomalous conservation equation 

d ,ja = ±e^F, u Va. (5) 

As in the electroweak theory, the anomaly appears in the divergence of a vector current due 
to the asymmetry in couplings between left- and right-handed fermions to the gauge field. 

3 CP violating phases will be ignored. 
An action with canonically normalized kinetic terms but rescaled parameters is obtained by rescaling 
the bosonic fields by 1/yF. 
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The vacuum structure of this theory is non-trivial and is described by the homotopy 
group, 7Ti(S' 1 ) = Z, as it is in the electroweak theory. Each inequivalent vacuum is labelled 
by the Chern-Simons number: 

N cs = I dxe 0u A» (6) 

which takes values in Z. Imposing periodic boundary conditions in a box of size L, the 
classical vacuum with topological number Ncs is given by 

* = *e*0»> , A, = ^d.aix), (7) 

where a(0) — a(L) = 2n N C s and N C s G Z. Working in the temporal gauge A = 0, the 
Coulomb gauge condition d\A\ = fixes the remaining time- independent gauge freedom and 
eq.| determines the vacuum gauge field to be {A — 0, A\ — 2ir N C s/eL}. 

Fermion number violation occurs for each fermion flavor when the spacetime integral of 
eq.^is non-vanishing. This is the case for gauge fields with non-vanishing winding number, 
defined as 

v = ^ J d 2 xe^F» v . (8) 

Configurations with non-zero winding number describe transitions between adjacent vacuua 
since 

N cs {t = +oo) - N cs (t = -oo) = ^- {dtdx — e 0v A v = v. (9) 

27T J at 

Transitions between adjacent vacuua, whether by tunnelling or by finite temperature fluctu- 
ations, are necessarily accompanied by fermion number violation, according to eq.|5] and eq.[| 
At zero temperature, the Abrikosov-Nielsen-Olesen vortices are Euclidean configurations, or 
instantons, which describe such tunnelling events JT^]. At finite temperature, sphalerons, to 



be discussed below, are assumed to be the relevant configurations for transitions between 
adjacent vacuua. 

In the absence of fermions, it is well-known that the 1+1 ABH model does not exist 
in the Higgs phase. Instantons destroy the long range order and a confining phase arises 

7 



17| . However, the existence of fermion zero modes in the instanton background implies the 



vanishing of instanton contributions to fermion-number-conserving Greens functions ||18||. In 
particular, the effective potential which determines the ground state of the theory receives 
no contribution from instantons, and there is no restoration of symmetry when fermions are 
included in the model|18|, [HJ. Therefore, this model is expected to exhibit the symmetry- 
breaking implied by the tree level potential. 

The 1+1 ABH theory contains a non-contractible loop in the configuration space of the 
bosonic sector of the theory which interpolates between vacuua differing by one unit of 



iV(7ff ||14|| . The field configuration of maximum energy on this loop is the sphaleron, which 
determines the height of the free energy barrier between adjacent vacuua in the broken 
phase of the theory. The gauge field on the non-contractible loop has the vacuum form 
{A = 0,A 1 = 2 7rN cs (r)/eL}. 

The scalar field on the non-contractible loop is now simply described by a scalar field 
in the form = e - 2mN cs(-r)x/L which removes the dependence on the constant 

gauge field|T^]. Identifying the Chern-Simons number on the loop, Ncs(t), with the loop 
parameter itself r = [0, 1], and solving the time-independent equations of motion then gives 
the scalar field : 



V \ iP' V 

Re& (x) = — ^=tanhy— ^— x, Jm$ (x) = — ^= cos it t, (10) 

where v 2 = v 2 sin 2 nr. By convention, the Imaginary part of $' has been chosen to be 
spatially constant. So, the Real part of $' interacts in a "slice" of the Higgs potential, which 
is a double-well potential. The solution is then a familiar "kink" configuration of the theory 
of a single scalar field interacting through a double-well potential, with a mass controlled by 
the loop parameter r. 

The energy functional reduces to E{r) = 2 -rp \^Xv 3 | sin 3 tct | for configurations on this 
loop. The energy reaches a maximum when r = 1/2, and this defines the unstable sphaleron 



configuration. In terms of the original variables, it takes the form: 



$ = ~J2 tanh \\^~T X ) ex P( i7ra; / L ) ' Ai = 7l j ^o = 0. (11) 

This is an unstable static solution to the source-free equations of motion for the scalar field. 
The r-direction of the energy functional along the loop represents the only unstable direction 
in configuration space away from the sphaleron. The sphaleron energy, E sp h = y/\v 3 = 
| v 2 M^, has the familiar form of a non-perturbative result since the weak coupling regime 
is achieved for v » 1 0. In the following section, the contribution of fermions to the free 
energy of this system will be calculated to determine thermal and quantum corrections to this 
classical sphaleron configuration. The constant gauge field in eq.[TT| does not contribute to 
the fermion determinant to be calculated, and so the gauge field will be neglected below [ PH . 



In this sense, the approach in this paper is complimentary to where the Yukawa coupling 
is neglected but bosonic fluctuations around the sphaleron are calculated. 

The sphaleron plays a key role in determining the rate of fermion number violation in the 
theory. As determined above, the sphaleron energy is the height of the free energy barrier 
between adjacent vacuua. If the sphaleron is also assumed to be the unique saddlepoint 
of lowest free energy, then its energy controls the rate of finite temperature fluctuations 
between adjacent vacuua. As mentioned in the Introduction, this rate is accompanied by a 
Boltzmann suppression factor 

r cxexp[-(3E sph ] (12) 

for temperatures much lower than the barrier height but much larger than the scalar mass, 
Mff, << T << v 2 M ( p. Similarly, the rate of fermion number violation which accompanies 
vacuum transitions is also determined by the Boltzmann exponential in this temperature 
range The strong dependence of the rate on the energy of the sphaleron requires an 
accurate determination of this quantity and provides further motivation for the work in this 
paper. 

5 Analogously, the electroweak sphaleron has energy, E sp h oc M w /a w . 



3 The Free Energy in the WKB Approximation 

The partition function for the system of fermions in a background scalar field, $, can be 
represented Euclidean path integral, 



Z [$] = exp [-N In det (i ft E - gv)] J iDtpDip^ exp - J dt J d 3 xC E + 5C 



(13) 



The fermion fields ip and ^ are required to be anti-periodic with period (3 = 1/T. Ce is the 
Wick rotation of the Lagrangian given in eq.|], and 5C is the counterterm Lagrangian. The 
determinant factor is the finite temperature determinant in the constant background which 
serves to normalize the free energy. The counterterm Lagrangian 5C is given by 

SC = AN(\$\ 2 -v 2 /2), (14) 

and is sufficient to renormalize both the one point and the two point function. A is chosen 
so that the tadpole graph vanishes when $ = v. This prescription gives 

g 2 f°° dp 

2^ Jo vV + g*v 2 ' (15) 
Using this value for A, the physical mass will be 

m 2 phys = N(-2\v 2 -g 2 /7r). (16) 

This counterterm is sufficient to renormalize the theory at finite temperatures. 

When the gaussian integral over the fermionic fields is performed, the partition function 
is expressed as the determinant of the operator 

en 

where ui n are the Matsubara frequencies given by ui n = 2n(2n+l)/ (3 and e are the eigenvalues 
of the spatial operator O given by 

O = 7o B [i - \<PP L - \*<j)*P R ] . (18) 
10 



Pl and Pr are the projection operators on the left and right moving states, respectively. 

The free energy of this system, defined by the usual thermodynamic relation, F = —-khiZ, 
is an effective action for the scalar field, <fi. The remainder of the paper will focus on a cal- 
culation of this quantity using the WKB approximation. 

Performing the sum over the Matsubara frequencies in eq.[T7] leaves f] 

F = "T EE [/(<0- /(<*)]■ ( 19 ) 
where /(e) is the contribution to the free energy from states with energy e 

/(e) = e + | In (l + e"*) . (20) 

The sum over a corresponds to the sum over two helicity states, and eo are the eigenvalues 
of the operator O with <fi taking on the constant value v. Eq.lHJ] reduces to the shift in the 



Dirac sea in the zero temperature limit p2 



The calculation of the free energy has been recast as a calculation of the spectrum of 
the operator O. In terms of the up (+) and down (-) components of the wave function, the 
eigenvalue problem may be expressed as a Schrodinger equation, 

<^±-( R i + ! ^f + d ' ai ) 1 , ± = -El i , ± . (21) 
ax z \ ax ax j 

A convenient basis for the two dimensional Dirac matrices has been chosen in which 70 = o\ , 

7i = £03 , and 75 = 02. The scalar field has been expressed here in terms of R = Re(g $) and 

I = Im(g$) Q. 

The imaginary part of the scalar field will be assumed to play no role in the energetics 
of the sphaleron. This is clearly seen in the temporal gauge where the imaginary piece is 
utilized only to enforce the boundary conditions at x — ±00. The gradient energy of the 
phase in eq.[ll] will vanish in the infinite volume limit. Thus, the eigenvalue equations eq.|21 

can be decoupled in the temporal gauge for the general sphaleron- like background. 

6 For notational convenience we don not include here the ideal gas contribution 
7 We will assume that there are no CP violating phases. 
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It is convenient to follow the notation used in ref. |12| and introduce the parameters 



Xd = \\v^ ' V = ' Z = X l Xcl ' ^ 

The Compton wavelength of the scalar, x c i, and the inverse temperature, /3, are now the 
only dimensionful parameters. The size of the classical sphaleron is roughly of order x c \. 
The parameter, y, is the ratio of the scalar and fermion Compton wavelengths. Now the 
eigenvalue equation may be written in the simple form: 

dz 2 



- y*V ± {z) - + xlel 



V± = 0, (23) 



with a "potential" given by 

V4W = W - 1 T £ (f z ) • (24) 

in terms of a canonically-normalized Real scalar field, = \/2 Re(&). 

The eigenvalue problem will now be solved using a WKB approximation [D], p^ |. This 
approximation is expected to be valid when the scalar field varies on scales larger than the 
fermion Compton wavelength, or roughly for large y. Since the goal is to study the effects 
of fermions on the sphaleron, only background scalar fields with the same topology as the 
sphaleron will be considered. These consist of fields which interpolate between two minima 
of the potential as x varies from — oo to +oo, as in eq.|TT|. For these fields, the potential V 
will be parity even, and the eigenstates will be classified by their parity quantum number. 
Furthermore, these eigenstates naturally fall into two classes: continuum states and discrete 
(bound) states. 

3.1 The Continuum States 

The asymptotic form of the WKB wavefunctions will be given by plane waves 

rj en (k,z) - cos(W^r n ) , rf d (k,z) - sin(^±^r) , *^oo. (25) 
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The momenta here are given by the free dispersion relation: 



k = \lx 2 cl e 2 — y 2 



(26) 



Periodic boundary conditions in a finite box of spatial size L require that the WKB 
phase shifts 5 a be related to the momenta by: k° an L + 5 a (k) = 2irn. As the size of the box 
is eventually to be taken infinite, eq.|l9| can be expanded in powers of S/L, leaving 

5 2 > 



-jEE 

L <T k 



d£de_ 
de dk 



L 



O 



k J \ Li ) \ L 2 j 

where the phase shifts for even and odd parity wavefunctions have been combined, 



(27) 



r / jreven , rodd 



(2? 



The phase shifts contain all of the information about the potential and are given by the 
WKB approximation. 



5 a (k) 



dz [k a (z) - k] , k„(z) = \jx 2 d t 2 -y 2 - y 2 V a (z). 



(29) 



In the limit of large L, the sum in eq.|27| becomes an integral and the contribution of 
continuum states to the free energy is: 

2 



CO) it 



2ttx c1 V ^WTy 1 



1 - 



(30) 



1 + e^( fc ) 

with e given in terms of k by the free particle dispersion relation, eq.p6| 

The first term in the square brackets in eq.^ is the zero temperature contribution. It is 
UV divergent as it stands and must be regulated, after which the momentum integral can be 
carried out explicitly. The counterterm discussed previously, eq.TE, renormalizes this term 
and leaves fl2f 



F° 

cont 



y 2 N 



Aixx 



cl 



E 



dz 



2 I v 2 . 



-V a 



l + V a ) In l + \ -V a 



(31) 
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The second term in the square brackets in eq.^ is the finite temperature contribution 
and may be written in terms of the ratio of temperature to fermion mass, a = (3y/x c i. 

F(a\ ont = - ^ £ r dz r % X L= { - a 2 ^) - ) (32) 

vr/? 2 V -/-oo •'0 v^ 2 + a 2 1 + e Vx +a I v J 

This integral is UV finite, but cannot be evaluated analytically for arbitrary values of the 
temperature-to-mass ratio, a. However, its limiting forms can be. The finite temperature 
contributions will be evaluated in the limit of large and small a? below. But first, the 
contribution of discrete states to the free energy must be included. 

3.2 The Discrete States 

The spectrum of discrete states is determined in the WKB approximation by the Bohr- 



Sommerfeld condition. The Schrodinger equation, eq.|21], will have a discrete bound state 
whenever 

w a (e) = - / dzk„(z)e(kl(zj) . (33) 



7T J-oo 

equals half of an odd integer. The total number of discrete states is then given by the integer 



closest to w a {y/x c i). Therefore, the contribution to the free energy, eq.[19|, from the discrete 
states will be given by (12] 

«r w a (y/x c i) 

F disc = --Y: E f(e)~ f(y/x cl ). (34) 

A cr to CT >0 

This sum may always be expressed in terms of an integral 

F dlsc = -- £ ^ dw [/(e) - f(y/x d )] + Rem. (35) 

with a remainder given by the Euler-MacLaurin formula. The remainder represents the error 
in the free energy from an integral approximation. As the "potential" becomes deeper, the 
number of bounds states w(y/x c [) increases, and the integral approximation becomes better. 
The depth of the potential is controlled by the dimensionless parameter y, the ratio of the 
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fermion and scalar Compton wavelengths. Thus, the integral approximation is obtained 
in the WKB regime, where the scalar field is slowly varying on the scales of the fermion 
Compton wavelength. This approximation will be discussed further below. 

Finally, eq.[35] may be cast in the same form as the contribution from the continuum 
states, eq.p2|, by changing variables from w to e. 



disc 



2tt V 



dz j 

-oo Jy^/l+V a /x c i 



de 



l + eP e 



X 2 i€ 2 _ y 2 {l + Va) 



(36) 



The first term in the square brackets is the zero temperature contribution. The integral 
over e can be done explicitly to give: 

v 2 n 



disc 



E 



dz 



-V a + 1 + K In 



(37) 



1 + V=K, 

The second term in the square brackets is the finite temperature contribution and may 
be written in terms a, defined previously. Changing variables to u = j3e gives: 



F(a 2 



'disc 



N 



7T 



E 



dz 



du 



aVT+v; 1 + e u 



a 2 (l + V a ) 



(38) 



Finally, the sum of the two zero temperature pieces, eq.|3T| and eq.BTl and the two finite 



temperature pieces, eq.|32] and eq.|38|, yields the total free energy, eq{L9 



3.3 The High Temperature Expansion 

The high temperature limit will provided a check on the previous expressions. One expects 
the leading finite temperature effects to be summarized in a temperature-dependent mass 
term for the scalar field. The high temperature expansion is defined as a 2 << 1. The free 
energy is not an analytic function of the temperature in this limit, and special care must be 
taken to obtain the expansion for small a (see Appendix). Including only the first non-trivial 
dependence on a 2 , eq.^ is: 
y 2 N 



F(a 2 ) cont 



2-RXr 



E 



dz 



- ~ (1 + V a {z) ) (2 7 - In (a 2 /n 2 ) ) 



1 + V a (z))ln 1 



0(a 2 



(39) 
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When this is combined with the zero temperature contribution, eq.^T], the square root and 
logarithmic terms cancel precisely. Then, the contribution of the continuom states in the 
high temperature limit is: 



F(a 2 ) cont = f^-E / dz 



1 r l-<j?/v 2 ) - l -V a {z) ( 2 1 -l + \n f(r> 



7T 2 



(40) 



The field gradients contained in V will cancel in the sum over a, leaving only a temperature- 
dependent mass term for the scalar field. Thus, it is clear that the finite temperature effects 
have made the theory local, in the sense that there is a well-defined gradient expansion. 



The contribution of discrete states follows from eq.|3q. It is not hard to demonstrate 
that the Euler-MacLaurin remainder is suppressed by at least 0(a). Furthermore, since the 
integral in eq.^ is cut off at the scale y/x c i, the Boltzmann exponential may be expanded 
in the high temperature approximation, and the zero temperature contribution cancels out, 
leaving 

F disc = 0(a). (41) 

As expected, low-lying bound states do not contribute to the free energy for high tempera- 
tures. 



Combining eq.[i(] and eq.|41] the total free energy in the high temperature approximation 
is given by 

We may now calculate the critical temperature in the high temperature approximation by 
combining this result with the tree-level scalar action from eq.[|. We will define critical 
temperature as the temperature for which the curvature at the origin vanishes. This is not 
necessarily the temperature at which the phase transition occurs. If the transition is first 
order, then bubble nucleation will induce the phase transition at a temperature above the 
critical temperature as defined here. Therefore, the critical temperature may overestimate 
the temperature at which the transition occurs, but this will be sufficient for this paper. 
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Requiring F tota i = at the origin gives 

(43) 

Interestingly enough, there is a range of temperature available in which a high temperature 
expansion, T >> y/x c i, is consistent with a broken phase, T < T c , and the sphaleron effects 
are of interest. So, for at least some range of heavy fermion mass, the leading temperature- 
dependent effect on the sphaleron is summarized by introducing a temperature-dependent 
mass for the scalar field, which in turn implies a temperature-dependent sphaleron mass. 
Thus, the sphaleron mass decreases with the scalar mass as the critical temperature is 
approached from below. However, notice that in the region of validity of the high T expansion 
(v » y), it is not expected that fermions will have a large effect on the sphaleron mass. 
Contrary to the hihg T case, in the low T limit the broken phase can be maintained for a 
wider range of heavy fermion masses. Therefore, we would expect that in this limit fermions 
will play a more important role in determining the sphaleron free energy. 

3.4 The Low Temperature Expansion 

The low temperature expansion is defined as a 2 » 1. It can be developed from eq.|32| by 
changing variables from x to u = (3t = \Jx 2 + a 2 . The finite temperature contribution to the 
free energy from the continuous states is then: 

^U, = - ^ E £ * f { + - } • (44) 

In the limit that the temperature is much less than the fermion mass, the Boltzmann dis- 
tribution is dominated by low energy states and the exponential in the integrand dominates. 

F(a 2 ) cont = - ^ E r dz due ~ U I \A 2 - a2 ( X + V °) - V^a 2 "} . (45) 

Tit) J —oo J a I J 



T, 



-7/4 y 

7r e " — exp 



2ttv 2 



17 



The second term in the curly brackets is independent of V and is simply expressed in 
terms of the MacDonald function, K\(a). 

poo ffVjf 

/ due~ u Vu 2 -a 2 = aK^a) = J— (1 + 0(1/ a)) e~\ a»l. (46) 

J a V 2 

The first term in the curly brackets cannot be evaluated in terms of simple functions. 
However, an integration by parts gives the leading temperature dependence immediately. 



roo , roo 

/ due- u Ju 2 -a 2 {l + V) = a^-V e~ a + / due' 11 



U 



^u 2 - a 2 (l + V) 



a 



V-V(l + 0(l/a)) e~ a , a»l. (47) 



So, all contributions of continuum states to the free energy are exponentially suppressed, 

0(e _a ), when the fermion mass greatly exceeds the temperature. 

The discrete states will give the dominant contribution to the free energy in the low 



temperature regime. This contribution is obtained from from eq.[37| and eq.|38|. The finite 
temperature part is: 

Nxri x-^ r°° , r a du 



^ = " 2 if ? L dZ L^V TT^ V" 2 - ^ + V.) . (48) 
For some point in the z-integration, the lower bound on the w-integral vanishes. Therefore, 



a simple expansion of the exponential as in eq.|45| will not be valid at this point. This is 
precisely the case where the background scalar goes through a zero. However, the large 
factor a » 1 allows —V to be quite close to unity and still allow the expansion. One 
requires a(l + V) » 1 which is —V « 1 — 1/a. If —V is not within 1/a << 1 of unity, 
then the low temperature expansion of the exponential is valid. 

If — V goes to unity at some point z Q , the ^-integral can be partitioned into three re- 
gions: z < z — b, z Q — b < z < z + b, and z > z Q + b. If b is greater than 0(1/ a), then the 
contribution of the first and third regions is exponentially damped since — V is not near 



unity. Then, the finite-T contribution to the free energy in the low temperature regime is 
approximately: 

Nxn ^ r zo+b , r a du 



Ha'W = "2 ^ E jT, * jf^ H -'HI + K.) + O (e-) . (49) 



The total free energy is obtained by adding the T-independent contributions, eq. |31] and 

eq.[37j. The remaining integral contains a region of size 2b in which —V approaches unity. 
A numerical computation of the free energy in this regime will be discussed in the next 



section. The qualitative behavior of eq|l9] is however clear: the finite-T contribution to the 



free energy decreases as the fermion mass increases in the low temperature regime. 

4 The Thermal Sphaleron 

As discussed in the previous sections, the thermal sphaleron is determined by minimizing 
the free energy while imposing the appropriate boundary conditions. In this section, this 
calculation will be performed in the high and low T limits. 

In the high T regime, the boundary conditions are temperature dependent due to the 
temperature dependence of the scalar expectation value. The scalar part of the free energy 
can be expressed in the form 



N r°° , 
F = — dz 

X c i J-oo 



(50) 



with a temperature-dependent vev defined as 



v\T) = ^-f^(7 + ln^). (51) 

Fermionic fluctuations are not expected to effect the form of the classical sphaleron which 
is given by | 



x 



<f>(x) = v tanh ( — ) . (52) 

v X'd > 



s This was demonstrated in the case of the soliton at zero temperature in ref. |12[ 

19 



So, only the effect on the size of the sphaleron relative to the size of the classical sphaleron 
x c i will be considered here. The phase in eq.[Tl| may be ignored for this purpose. It is 
easily seen from eq.BUI that the thermal sphaleron, in the high T case, is given by a simple 
modification of the classical solution. One replaces the vev v in eq.[52] by v(T) determined 



by eq.|51|, including its appearance in the scalar Compton wavelength x c [. The free energy of 
the ansatz in the high T approximation is then given by 

F sph = x i^T (53) 

where a temperature-dependent scalar Compton wavelength has been defined for convenience 



x{T) = ixm ■ (54 > 

This is just the classical sphaleron energy with the parameters appropriately changed to 
include the finite temperature effects. The free energy of the sphaleron is an increasing 
function of the fermion mass. Furthermore, the free energy increases only logarithmically 
with the temperature in two dimensions. As expected, when v » y, the thermal sphaleron 
reduces to the classical sphaleron because in this regime the classical contribution from the 
scalar sector overwhelms the fermionic fluctuations. Unfortunately, as was shown in the 
previous section, the high T approximation will break down in the regime where quantum 
fluctuations become important (i.e. y > v). 

It is interesting at this point to compare the free energy of the thermal sphaleron to the 
free energy of the classical sphaleron including the one-loop fermion contribution, which is 
given by 

3y 2 



_ 4AV 

1 ' cl — T, 



1 + 



Airv 2 



7 + In - J 

7T / 



(55) 



This free energy is indeed slightly higher than the thermal sphaleron free energy, eq.|53 



as 



expected. Of course, this comparison should only be made at temperatures low enough that 
the temperature-dependent vev of the scalar field is still near v. 
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In general , it may or may not be simple to solve for the thermal sphaleron in the 
high T limit. This depends upon whether or not the corrections to the gradient terms 
are suppressed by finite temperature effects. In the Standard Model at one loop order, 
the gradient corrections are suppressed because there is no renormalization of the kinetic 



term This is also the case in the 1+1 Abelian Higgs model. However, in this model, 
the suppression persists to all orders in the loop expansion because only the mass and the 
vacuum energy require renormalization in two dimensions. The suppression of the corrections 
to the gradient terms allows one to determine the free energy functional by calculating the 
fermionic fluctuations in a constant background. The effect of the fermions is to give the 
couplings in the scalar sector temperature dependence. The thermal sphaleron is given by the 
classical sphaleron with its parameters appropriately replaced with temperature-dependent 
parameters. This is the standard procedure when working in the high T limit. f22[ 

In the low T limit, one is no longer bound to the regime where v » y and the effect 
of fermion fluctuations may become more important. In this regime, it is not possible to 
get a closed form expression for the size of the sphaleron. The thermal sphaleron will be 
determined by using the following ansatz for its shape 

z x 
<Pz (x) = v tanh( — ) , z = — , (56) 
Zq x c i 

and minimizing the free energy with respect to the parameter, zq. The use of this ansatz 
assumes that the temperature is small enough that the scalar expectation value is still near 
v. The free energy of this ansatz in the low T limit is given by 

y 2 N 



2V 2 ( 1 \ ( 3 7T 

Vr + -J +Z °{^- 8 



+ -Ar / dz ( ^ U \/u 2 — a 2 t&nh z (z / z ) 

7TGT J-oo Ja|tanh(z/zo)| e u + 1 



(57) 



The first term in this equation is the classical contribution, while the second term is the 
zero temperature contribution from quantum corrections. In the T = limit, the size of the 
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sphaleron which minimizes eq.|57] is [12 



3y 2 ( 3 7i 2v 2 
2V* \ 2^ ~ 8 + V 



) 



1-1/2 



(58) 



Note that at zero T ,for the choice v/y = .1, the sphaleron size changes by a factor of 1/4 
relative to its classical value. Thus the effect of quantum fluctuations alone will tend to 
increase the sphaleron mass as the fermion is made heavier. Furthermore, the temperature- 
dependent effects induced by the fermions are expected to become less important as the 
fermion mass increases, because thermal fluctuations will be suppressed. 

As in the high T case, we expect the classical contribution to dominate the contribution 
from fermionic fluctuations in the limit that v » y and the size of the sphaleron will reduce 
to the classical value x c i, or Zq = 1. However, allowing y to become large while holding the 
temperature fixed, fermionic fluctuations become increasingly important, and there are non- 
negligible temperature-dependent corrections. The free energy was calculated numerically 
for several values of Zq and a. The thermal sphaleron size was determined by minimizing the 
free energy with respect to Zq, for a given value of a. Figure 1 shows the thermal sphaleron 
size in units of the T = sphaleron mass, as a function of the parameter a for v/y = 0.1. 
This figure shows that in addition to the quantum fluctuations thermal fluctuations are 
non-negligible (yet smaller than the quantum effects) as the temperature is increased. Thus, 
fermionic fluctuations will tend to increase the free energy of the sphaleron relative to its 
classical value. Figure 2 shows the free energy of the thermal sphaleron in units of the 
classical sphaleron energy for several values of the parameter a, with the choice v/y = 0.1. 
Figure 2 indicates that quantum corrections may change the sphaleron free energy by a factor 
of 5. Furthermore, within the validity of the low T expansion, the temperature-dependent 
effects are at the level of 20 percent. 
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5 Discussion 



The 1 + 1 Abelian Higgs model has illustrated the effects of heavy fermions on a non- 
perturbative scalar background. In both the low and high temperature regimes, a self- 
consistent calculation of the free energy of the sphaleron has shown that heavy fermions will 
increase the mass of the thermal sphaleron. It has been shown that the thermal sphaleron 
has a smaller free energy than the sum of the classical sphaleron energy and its fermionic 
corrections. That is to say the classical sphaleron is not a true saddlepoint of the free energy 
functional. 

Beyond the implications this result might have for phenomenology, the results obtained 
here also address the interesting theoretical questions mentioned in the Introduction. In 
the low temperature regime, the free energy was found to be a non-local functional of the 
scalar field. The free energy is correctly expressed in this regime only by the summation 
of gradients to all orders, a summation provided by the WKB approximation. However, 
the net effect of the fluctuation-induced gradients is suppressed exactly when one would 
expect quantum fluctuations become important (y > v). An heuristic explanation is that 
the sphaleron size acts as a cut-off for infra-red divergences. Since the sphaleron gets smaller 
when the fermionic mass (y) increases, the infrared cut-off increases, and fluctuation-induced 
gradients will be further suppressed. Of course, as the size of the sphaleron decreases, the 
kinetic term in the classical Lagrangian will become more important. Though the effects of 
gradients are small, thermal effects can still change the sphaleron free energy substantially. 
This is a consequence of the occupation of a large number of low-lying fermionic states 
which become available when the Yukawa coupling is increased. The T = quantum effects 
can also change the sphaleron free energy substantially as a consequence of the standard 
perturbative non-decoupling when the fermion mass is taken to infinity while holding the 
scalar vacuum expectation value constant. 

In the high temperature regime on the other hand, it has been shown here that the 
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gradient expansion is well-behaved. So, the net effect of thermal fluctuations at high T 
effects is the localization of the theory. As expected, the fermions do not decouple in this 
regime, even though their masses are determined by the dimensionful parameter, T, since 
they contribute to temperature-dependent renormalizations of the parameters in the theory 
which are physically observable. 

Thus, one may conclude that, in the simple model studied in this paper, the fermions 
do not decouple in the perturbative sense and can enhance the sphaleron mass. Since the 
gradients are suppressed, one may also conclude that the effects of the non-perturbative 
variation of the field discussed in the Introduction are small. One would expect similiar 
behaviour also in four dimensions for the physical reasons discussed above. The fact that 
there is only perturbative non-decoupling also encourages one to expect that non-perturbative 
methods may not be necessary in order to calculate the effects of heavy fermions on non- 
perturbative scalar backgrounds in four dimensions. 

6 Appendix 

The integral, eqj32], is not an analytic function of a in a neighborhood of a = 0. So, a 
naive Taylor expansion in the region, a ~ 0, is bound to fail. This is a familiar problem 
which occurs in the calculation of the free energy of a gas of free particles. The work in this 
Appendix will extract the leading non-analytic behavior near a ~ 0. 
Consider the integral: 



where x = kf3/x c i and a 2 = (3 2 m 2 = (3 2 y 2 v 2 / x 2 cl . The integral is finite at a = 0, reproducing 
the free energy of a gas of free fermions at high temperatures: 




(59) 



I(0;V) = 2-T 2 . 



(60) 



This high temperature result is independent of the background field, as expected. 
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The first derivative of eq.|59| is: 
dl , 2 TrN v 7 ^ 7 1 , y /-\ f°° dx 2 



da^ a]V) ea + l 2 (1 + *° Jo ^2 _ fl2y) (x2 + fl2) + x ( 61 ) 

after an integration by parts. The surface term is 0(a°) but the integral is logarithmically 
divergent as a 2 — > 0, so one can expect: 

dl 1 

— ( a 2 ; V ) V^V + const, log a (62) 

oa z 2 

near a 2 = 0. 

To determine this non-analytic term, the integral is broken into a series by use of the 
identity: 

1 1 +oc y 

ey + 1 ~ 2 ~ ^ 2/ 2 + ?r 2 (2n + l) 2 ^ 

which follows from a contour integration of the left hand side [[2^]. Now, the integrals over 

each term in this expansion are divergent, so they will be regulated by introducing: 

/>oo A ™— € 

Je(a 2 ;V)= . (64) 

JO ^ x 2 _ a 2y^ ( x 2 + fl 2) e^+^ + 1 



and 



where 



J £ (a 2 ; V) = J £ (1) (a 2 ; V) + J £ (2) (a 2 ; V). (65) 



and 



J« (a 2 ; V0 = - V / dx 2 . :- ^ ^ - (66) 

e V ; ^ Jo v / x 2 -a 2 V r x 2 + a 2 + 7T 2 (2n+l) 2 V ; 

j(2) (a 2 y) = - /°° dx 2 X ~ 6 . (67) 

A little analysis shows that and J^ 2 - 1 converge if < e < 2. The integrals will be 
estimated for small a 2 in this region, and then analytically continued to e = 0, which is a 



regular point of the original integral, eq.64. The result is the unique high-T expansion of /. 
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First, eq.^6] can be rewritten as: 



+00 

E 



x 



» [n 2 (2n + l) 2 ] 1 ^ [1 + a 2 /7r 2 (2n + l) 2 ]^ 

du = 

(M 2 + l) v /u 2 + c 2 (n) 



(68) 



where c 2 (n) is an n-dependent constant of order a 2 . To leading order in a 2 then, eq. 38 is: 

1 



jP(a 2 ;V) = -£ 



l + e 



1 + 0(a 2 



(69) 



^ [vr 2 (2n+ I) 2 ]"* - J o u{u 2 + l) 

The remaining integral is elementary, and the sum can be expressed in terms of the Riemann- 
Zeta function. 



JP (a 2 ; V) = - 2 (l - 1/2^) C [1 + e] £^ + 0(a 



(70) 



Finally, this result is continued to e = by means of the formula 

2 e 7r 1+e C[-e] 



C[l + e] 



sin(f ) r [1 + e] 



which gives 



jW ( a 2 ;V) = -1/e - 7 + log(7r/2) + 0(e) + 0(a 2 ) 
where 7 ~ 0.53 is the Euler Number. 



(71) 



(72) 



Next, eq.57 can be rewritten as: 



JW (a 2 ; V) = a-*jP (1; V) = ( 1 - eloga + 0(e 2 ) ) (1; V) 



(73) 



This is the source of the non-analyticity near a = 

Now, J™ (1; V) is easily analyzed by expressing it in terms of two special functionspO 

1 



4 2) (l; V) 



2V-F 



5(1 - e/2, e/2) 2 F X [1/2, 1 - e/2; 1; 1 + l/V] 



(74) 
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The pole at e = is contained in the Beta function, B(l — e/2, e/2). The hypergeometric 
function, 2 -Pi, is analytic in its second argument near 1. Expanding in a Taylor series there 
gives 



JP(1;V) = i 5(1-6/2, 6/2) 



1 



2 F 1 '[1/2,1;1;1+1/V] + 0(e 2 



(75) 



2 V /Z F 

where the prime denotes differentiation with respect to the second argument of 2-^1- This 
function can be determined by using the various recursion relations which relate hypergeo- 
metric functions of different arguments. The result is: 



iF[ [1/2, 1; 1; 1 + 1/V] = 2 V^V log 



- (1 + V^v 



Now, expanding the Beta function for small e in eq.[74] yields: 



J®(1;V) = 1/6 - log[- (l + 



+ 0(6). 



Finally, eq. 73 becomes 



J® (a 2 ; V) = 1/6 - log(a) - log [- (l + 



+ 0(6) 



Combining eq.ff^, eq.|78] and eq.|65], the poles in e cancel and give 



J e (a 2 ; V) = - log(2a/7r) - j E - lo. 
dl 



+ 0(a) + 0(e) 



(76) 



(77) 



(7S 



(79) 



(a 2 ;V) 



\ V^V + I (1 + y) I log(2a/7r) + lE 



0(a) + 0(e) . (80) 



da? x ~ 1 ' ' 2 

+ log [1 (l + 

Now, /(a 2 ) can be recovered by expanding the surface term in eqj61| to lowest order in a, 
combining with eq.[8Fj, and integrating with respect to a 2 . 

I(a 2 ;V) = I(0-V) l a 2 V-V+ \ (l + V) a 2 { log(a 2 /7r 2 ) + 2 7s - 1 } 

+ 0(a 4 ) (81) 



2 
1 



+ - (l + V) a 2 log[l + V^V 

where the integration constant is determined by eq.|K]. This result establishes the free energy 
in eq.|39. 
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Figure Captions 

Figure 1: The ratio of the size of the thermal sphaleron to the zero temperature sphaleron 
as a function of the parameter a = j3y/x c i. 

Figure 2: The free enrgy of the thermal sphaleron as a function of the parameter a = /3y/x c [. 
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